Infrared behavior of interacting bosons at zero temperature 
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We review the infrared behavior of interacting bosons at zero temperature. After a brief discus- 
sion of the BogoUubov approximation and the breakdown of perturbation theory due to infrared 
divergences, we present two approaches that are free of infrared divergences - Popov's hydrodynamic 
theory and the non-perturbative renormalization group - and allow us to obtain the exact infrared 
behavior of the correlation functions. We also point out the connection between the infrared behav- 
ior in the superfluid phase and the critical behavior at the superfluid-Mott-insulator transition in 
the Bose-Hubbard model. 
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Many of the predictions of the Bogoliubov theory of 
superfluidity [1] have been confirmed experimentally, in 
particular in ultracold atomic gases Nevertheless a 

clear understanding of the infrared behavior of interact- 
ing bosons at zero temperature has remained a challeng- 
ing theoretical issue for a long time. Early attempts to go 
beyond the Bogoliubov theory have revealed a singular 
perturbation theory plagued by infrared divergences due 
to the presence of the Bose-Einstein condensate and the 
Goldstone mode 0-0|- 1970s, Nepomnyashchii and 

Nepomnyashchii proved that the anomalous self-energy 
vanishes at zero frequency and momentum in dimension 

^ 3 This exact result shows that the Bogoliubov 
approximation, where the linear spectrum and the super- 
fluidity rely on a finite value of the anomalous self-energy, 
breaks down at low energy. As realized latter on [ld| , 
the singular perturbation theory is a direct consequence 
of the coupling between transverse and longitudinal fluc- 
tuations and reflects the divergence of the longitudinal 
susceptibility - a general phenomenon in systems with a 
continuous broken symmetry 

In this paper, we review the infrared behavior of inter- 
acting bosons. A more detailed discussion together with 
a comparison to the classical O(A^) model can be found 
in Ref. [l^. In Sec. m we briefly review the Bogoli- 
ubov theory and the appearance of infrared divergences 
in perturbation theory. We introduce the Ginzburg mo- 
mentum scale PG signaling the breakdown of the Bogoli- 
ubov approximation. In Sec. IIII[ we discuss Popov's hy- 
drodynamic approach based on a phase-density repre- 
sentation of the boson held flMl. This approach al- 
lows one to derive the order parameter correlation func- 
tion without encountering infrared divergences. The non- 
perturbative renormalization group (NPRG) provides an- 
other approach free of infrared divergences which yields 
the exact infrared behavior of the normal and anomalous 
single-particle propagators fSec. lIVp (l6l-[2^. In the last 
section, we discuss the critical behavior at the superfluid- 
Mott-insulator transition in the Bose-Hubbard model 
and its connection to the infrared behavior in the su- 



II. PERTURBATION THEORY AND 
BREAKDOWN OF THE BOGOLIUBOV 
APPROXIMATION 

We consider interacting bosons at zero temperature 
with the (Euclidean) action 



S 



dx 



2m / 2 



(1) 



where ^{x) is a bosonic (complex) fleld, x — (r,T), and 

J dx — Jq dr J d'^r. r S [0, /3] is an imaginary time, (3 — > 
oo the inverse temperature, and denotes the chemical 
potential. The interaction is assumed to be local in space 
and the model is regularized by a momentum cutoff A. 
We consider a space dimension d > 1. 
Introducing the two-component field 



V'b) 

n-p) 



vl/t(p) = (V/(p),V(-p)) (2) 



(with p — (p, iLu) and a; a Matsubara frequency), the one- 
particle (connected) propagator becomes a 2 x 2 matrix 
whose inverse in Fourier space is given by 
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-San(p) 
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where En and San are the normal and anomalous self- 
energies, respectively, and Cp — p^/2m. If we choose 



the order parameter {tpiyx)) — y/riQ to be real (with tiq 
the condensate density), then the anomalous self-energy 
Sanb) is real. 



A. Bogoliubov approximation 

The Bogoliubov approximation is a Gaussian fluctu- 
ation theory about the saddle point solution ip{x) = 



2 





FIG. 1: One-loop correction S'^^ to the self-energy. The dots 
represent the bare interaction, the zigzag lines the order pa- 
rameter ^/no, and the solid lines the connected propagator 



/no = y^fijg. It is equivalent to a zero-loop calcula- 
tion of the self-energies |25l [26| , 



This yields the (connected) propagators 



(4) 



gno 



El 
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where Ep = [ep(ep + 2(/no)] ' is the Bogoliubov quasi- 
particle excitation energy. When |p| is larger than the 
healing momentum pc — (2(?m7io)^^^, the spectrum Ep ~ 

+ gno is particle-like, whereas it becomes sound-like 
for IpI <C Pc with a velocity c = ^JgriQ/m. In the weak- 
coupling limit, rto — n. (n is the mean boson density) and 
Pc can equivalently be defined as pc — (2gmny^^. 

It is convenient to write the boson field 



1 



V2 



[-01(2;) -I- i'ijj2{x)] 



(6) 



in terms of two real fields -01 and which allows us 
to distinguish between longitudinal (^i) and transverse 
(^2) fluctuations. In the hydrodynamic regime |p| <^ Pc, 



g'^,\p) = 



2gno 



(7) 



where Gy (p) — ('i/'i(p)V'j (~P))c- In the Bogoliubov ap- 
proximation, the occurrence of a linear spectrum at low 
energy (which implies superfluidity according to Lan- 
dau's criterion) is due to Ean(p = 0) being nonzero. 



B. Infrared divergences and the Ginzburg scale 

Let us now consider the lowest-order (one-loop) cor- 
rection E^^) to the Bogoliubov result S*^°'. For d < 3, 
the second diagram of Fig. [T] gives a divergent contribu- 
tion when the two internal lines correspond to transverse 



fluctuations, which indicates a breakdown of perturba- 
tion theory and therefore the Bogoliubov approximation. 
Retaining only the divergent part, we obtain 

E(i)(p) ^ S«(p) ^ -iff^no^G(°Hg)G(?(p + g), (8) 
where we use the notation f — f°° ^ f . For small 

Jq J— 00 Ztt Jq 

p, the main contribution to the q integral in ([8]) comes 
from momenta |q| < pc and frequencies |w'| < cpc- We 
can then use ([7]) and obtain 



eW(p)^eW(p): 

if (i < 3 and 

eW(p)^eW(p) 

if c? = 3, where 



-2^A,+r ( p^ 



2\ (^-3)72 
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p2 + UJ^/c^ 



-<i/2 r(d/2) 



Ad=< sin(7rd/2) r(d-l) 

I A if d 



if d < 4, 
4. 



(9) 
(10) 

(11) 



We can estimate the characteristic (Ginzburg) momen- 
tum scale pg below which the Bogoliubov approximation 
breaks down from the condition |si^''(p)| ^ I^n \p) or 

|Si'nHp)| 



Si°^(p) for IpI = Pg and |a;| = cpG, 



PG 



(^<j+i.gmp,)i/(3-rf) if d < 3, 
[-^r^] if d-3. 



Pc exp 



A4gmpa ^ 

This result can be rewritten as 

p,(Arf+ig'^/2)i/(3-<i) if < 3, 

P^^^P(- A, 7^33/2 ) if d = 3, 



PG 



where 



-1-2/rf ^ 



(12) 



(13) 



(14) 



is the dimensionless coupling constant obtained by com- 
paring the mean interaction energy per particle gn to the 
typical kinetic energy 1 /mf^ where f ~ n-i/'' is the mean 
distance between particles [23|. A superfluid is weakly 
correlated if 5 ^ 1, i.e. pg ^ Pc ^ n}^'^ (the character- 
istic momentum scale n^l'^ does however not play any role 
in the weak-coupling limit) [2^. In this case, the Bogoli- 
ubov theory applies to a large part of the spectrum where 
the dispersion is linear (i.e. |p| < p^ and breaks down 
only at very small momenta |p| ^ pg ^ Pc- When the di- 
mensionless coupling g becomes of order unity, the three 
characteristic momentum scales pG ^ Pc ^ n^l'^ become 
of the same order. The momentum range \pgiPc\ where 
the linear spectrum can be described by the Bogoliubov 
theory is then suppressed. We expect the strong-coupling 
regime g 3> 1 to be governed by a single characteristic 
momentum scale, namely n^/''. 
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C. Vanishing of the anomalous self-energy 

The exact values of Sn(p — 0) and San(p — 0) can be 
obtained using the U(l) symmetry of the action, i.e. the 
invariance under the field transformation ^/'(cc) — >■ e^^ip{x) 



and iIj*{x) 



-ie 



i^*{x). On the one hand, the self- 



energies satisfy the Hugenholtz-Pines theorem [6|, 

En(p = 0) - Ea„(p = 0) = (15) 

On the other hand, the anomalous self-energy vanishes, 



Sa„(p - 0) = 0. 



(16) 



The last result was first proven by Nepomnyashchii and 
Nepomnyashchii fl,[H,[2J|- It shows that the Bogoliubov 
theory, where the linear spectrum and the superfluidity 
rely on a finite value of the anomalous self-energy, breaks 
down at low energy in agreement with the conclusions 
drawn from perturbation theory (Sec. Ill Bp . 



where Ep is the Bogoliubov excitation energy defined 
in Sec. Ill Al In the hydrodynamic regime |p| ^ Pc = 



Gnnip) 



n p 



GnB (P) — J— : 



(20) 



mc 1 



where c = y^gh/m is the Bogoliubov sound mode ve- 
locity. It can be shown that Eqs. ((20|) are exact in the 
hydrodynamic limit |p|, |w|/c <C Pc provided that c is the 
exact sound mode velocity and n the actual mean density 
(which may differ from fi/g) [l^ [l3|. 



B. Normal and anomalous propagators 



III. POPOV'S HYDRODYNAMIC THEORY 

It was realized by Popov that the phase-density repre- 
sentation of the boson field i b = v/ne'^ leads to a theory 
free of infrared divergences |13l Il4| . In this section, we 
show how this allows us to obtain the infrared behavior of 
the propagators G„(p) and Gan(p) without encountering 
infrared divergences (isj . 



A. Hydrodynamic action 

In terms of the density and phase fields, the action 
reads 



S[n,e] 



dx 



Tl 

in9+—(V(, 
2m 



(Vn) 



, 9 1 



8mri 

(17) 

At the saddle-point level, n(x) = n = fi/g. Expanding 
the action to second order in Sn = n — fi, and V0, we 
obtain 



S[dn,9] = / dx 



iSn9 



(18) 

The higher-order terms can be taken into account within 
perturbation theory and only lead to finite corrections of 
the coefficients of the hydrodynamic action ([TSl) [14]. 

We deduce the correlation functions of the hydrody- 
namic variables, 



Gnnip) = {Sn{p)Sn{-p)) = 
Gneip) = {5n{p)e{-p)) = - 



£^2 
p 



(19) 



Geeip) = {e{p)e[-p)) 



^— + a 

Amn ^ 



To compute the propagator of the ip field, we write 

i^{x) = v/no + M^)e''^"\ (21) 



where hq = |(V'(a;))p = |(-\/n(a;)e'^^'^^)p is the conden- 
sate density. For a weakly-interacting superfluid, hq ~ h, 
and we expect the fiuctuations Sn to be small. Let us as- 
sume that the superfiuid order parameter {ip{x)) = y'no 
is real. Transverse and longitudinal fluctuations are then 
expressed as 



(5V2 = V^^o + 



2 



(22) 



where the ellipses stand for subleading contributions to 
the low-energy behavior of the correlation functions. For 
the transverse propagator, we obtain 



G22(p) ^ 2noGee{.p) = 



2nQmc^ 



1 



c2p2 



n oj^ 

to leading order in the hydrodynamic regime, while 
Guip) ^ GMp) = 



(23) 



(24) 



The longitudinal propagator is given by 



1 



no 



Giiix) = ^Gnn{x) + '^{e{xfe{Qf) 

Zhq Z 



2nQ 



Gnn{x) + noGgg{x)^, 



(25) 



where the second line is obtained using Wick's theorem. 
In Fourier space. 



Gii(p) 



2mno io^ -f c^p^ 



noGee*Gee{p), (26) 
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where 



Gee * Geeip) = / Gee{q)Gee{p + q). (27) 

The dominant contribution to the integral in (j27p conies 
from momenta |q| "^Pc and frequencies |w'|/c <Pc, i-e. 



Gee * Gee{p) 



^<i+ic(f) (p2 + ^) ifd<3, 



2 V fi , 



p2 + ^ 



if d = 3. 



(28) 



By comparing the two terms in the rhs of (j26p with 
IpI — Pg and |a;| = cpc, we recover the Ginzburg scale 
P^ . For IpI, |w|/c ^ the last term in the rhs of 
(PS)) can be neglected and we reproduce the result of 
the Bogoliubov theory (noting that n ~ tiq), while for 

IpI, |c^|/c < PG, Gii(p) ~ 1/(C^2 ^ ^2p2)(3-d)/2 ig doj^. 

inated by phase fluctuations (Goldstone regime). The 
longitudinal susceptibility Gii(p,icij = 0) ~ l/lpl^"*^ 
for p ^ in contrast to the Bogoliubov approximation 
Gii(p,itJ = 0) = 1/2toc2. 

From these results, we deduce the hydrodynamic be- 
havior of the normal propagator [2^, 



Gn(p) 



1 



[Gn(p)-2iGi2(p)+G22(p)] 



ngTOC 



n cj^ + c^p^ 



VjJ 



up' + c^p2 



:Gii(p), 



(29) 



as well as that of the anomalous propagator, 
1 



Gan(p) 



[Gn(p)-G22(p)] 



Home 



1 



n uj'^ + c^p^ 



Igii{p), (30) 



where Gii(p) is given by (|26p . The leading-order terms 
in and (150)1 agree with the results of Gavoret and 
Nozieres [TJ and are exact (see next section). The con- 
tribution of the diverging longitudinal correlation func- 
tion was first identified by Nepomnyashchii and Nepom- 
nyashchii ^] and later in Refs. [H-lB, [13, [3l| . 



C. Normal and anomalous self-energies 

To compute the self-energies Sn(p) and I]an(p), we use 
the relations 



S„(p)=Goi(p) 



San(p) 



Gn{~p) 



G„(p)G„(-p)-Ga„(p)2' 



Gan(p) 



(31) 



G„(p)G„(-p)-Ga„(p)2 



This yields QlTal 

San(p) = Sn(p)-Gol(p) 



j(cj2^cV)^^~''^/' if d<3. 



A4cnQm'^ 



In 



if d = 3, 
(32) 



in the infrared limit |p|, |w|/c <^ pq, where Gq^{p) = 
iu ~ ep + /i. Equations ([5^ agree with the exact results 



(|15ll6p and show that Sn(p) and San(p) are dominated 
by non-analytic terms for p — > 0. This non-analyticity 
reflects the singular behavior of the longitudinal correla- 
tion function 



Gii(p) 
in the low-energy limit. 



1 



2San(p) 



(33) 



IV. THE NON-PERTURBATIVE RG 

The NPRG provides another way to circumvent the dif- 
ficulties of perturbation theory and derive the correlation 
functions in the low-energy limit (l6l - [2^ . The strategy 
of the NPRG is to build a family of theories indexed by 
a momentum scale k such that fluctuations are smoothly 
taken into account as k is lowered from the microscopic 
scale A down to (3^ . [33| . This is achieved by adding to 
the action ([!]) an infrared regulator term 



(34) 



where ipi and ^p2 are the two real fields introduced in 
Sec. Ill Al [Eq. The main quantity of interest is the 

so-called average effective action 

rk[4>] = -In Zk[j]+Y,M--P)MP) - A5fc[0], (35) 



defined as a modified Legendre transform of — In Zk [J] 
which includes the subtraction of AS'fcIc/)]. Ji denotes an 
external source that couples linearly to the boson field -0^ 
and (f>i{x) = (ipi^x)) is the superfiuid order parameter. 
The cutoff function Rk is chosen such that at the micro- 
scopic scale A it suppresses all fluctuations, so that the 
mean-field approximation FaI^] — S[4>\ becomes exact. 
The effective action of the original model ([T]) is given by 
Ffe=o provided that Rk=o vanishes. For a generic value of 
fc, the cutoff function Rk{p) suppresses fluctuations with 
momentum |p| < k and frequency |w| < ck but leaves 
those with |p|, |cij|/c > k unaffected (c = Ck is the velocity 
of the Goldstone mode). The dependence of the average 
effective action on k is given by Wetterich's equation [33] 



dtVki^] = -Tr 



R, 



Ri 



(36) 
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where t = ln(fc/A) and Rk = dtRk- '[(j)] denotes the 
second-order functional derivative of rfc[0]. In Fourier 
space, the trace involves a sum over momenta and fre- 
quencies as well as the internal index of the ip field. We 
choose the cutoff function [2Q] 



Rk{p) 



2m 



£1 

fc2 



(37) 



where r{Y) = (e^ — 1)^^. The fc-dependent variable 
ZA.k is defined below. A natural choice for the velocity 
Co would be the actual (fc-dependent) velocity c^. of the 
Goldstone mode. In the weak coupling limit, however, 
the Goldstone mode velocity renormalizes only weakly 
and is well approximated by the fc-independent value 
Co = yjgn/m. 



A. Derivative expansion and infrared behavior 



Because of the regulator term ASfc, the vertices 
(-Pi;''' ,Pn) ^rc smooth functions of momenta 
and frequencies and can be expanded in powers of p| /fc^ 
and w^/c^fc^. Thus if we are interested only in the low- 
energy properties, we can use a derivative expansion of 
the average effective action (32l . [ssj . In the following we 
consider the Ansatz 



= I dx 



4>* (Zckdr - VamBI 



ZA,k 

2m 



(38) 



where n — \(j>\'^ — ^{(t>i+4i2)- ?^o,fc denotes the condensate 
density in the equilibrium state. We have introduced a 
second-order time derivative term. Although not present 
in the initial average effective action Fa , we shall see that 
this term plays a crucial role when d < 3 fisl, l2l'| . 

In a broken U(l) symmetry state with order parameter 
= \/2no, (/)2 = 0, the two-point vertex is given by 



^'k.LiP) = Vam^'' + ZA,kep, 



k,22 
.(2) 
fc,12 



(39) 



(p) = Zc,kUJ. 



Since these expressions are obtained from a derivative 
expansion of the average effective action, they are valid 
only in the limit |p|,|aj|/c<Cfc. In practice however, one 
can retrieve the p dependence of F^^q(p) at finite p by 
stopping the RG fiow at fc ^ + liP- jc? [20]. From 
pop we deduce the derivative expansion of the normal 
and anomalous self-energies, 

Sfc,n(p) = M + ^A.k^"^ + (1 - Zc,k)i'^ 

- (1 - ^A,fc)ep + >ykno^k, 

^k,an{p) = AfeTlo^fc. 



At the initial stage of the flow, Za.a = Zc.a = 1, Va,a = 
^1 — 9 and riQ^A = m/Si which reproduces the results 
of the Bogoliubov approximation. 

Since the anomalous self-energy 'S,k=a,a,n{p) ^ ("^^ + 
^2p2-j(3-ci)/2 jg singular for |p|, |w|/c ^ pa and d < 3 
(see Sec. mil) , we expect T,k,an{p = 0) ^ fc'^"'' for k <^pg 
(given the equivalence between k and -^/p^+w^/c?) , i.e. 

Afc - e-'' (40) 

(rto,fc=o is finite in the superfluid phase). The hypothesis 
(Uni) is sufficient, when combined to Galilean and gauge 
invariances, to obtain the exact infrared behavior of the 
normal and anomalous propagators. In the domain of 
validity of the derivative expansion, |pp, jw^/c^ <C fc^ <C 
^ q_ Qj^g obtains from ([59)) 



Gk,ii{p) 
Gk.,22ip) 

GkMp) 



1 



2Afeno,fc 
1 1 

VA.k + c2p2 ' 

Zc,k 



(41) 



2Afcno,fcVA,fe uj^ 



where 



Ck = 



ZAMf^m 
VA,k + Z^i,/2Xkno^k 



C^p2 



1/2 



(42) 



is the velocity of the Goldstone mode. From (HO)) and 
(|4T|) . we recover the divergence of the longitudinal sus- 
ceptibility if we identify fc with y/p"^ + w^/c^. 

The parameters ZA^k, Zc,k and YA,k can be related to 
thermodynamic quantities using Ward identities 0, ITtI . 

MM. 



na.k = ZA,kno.k = nk 
1 d^Uk 



VA,k 
Zc,k 



2^0, fe d^i'^ 
d^Uk 



(43) 



dndiJ, 



d^ ' 



where fik is the mean boson density and n^^k the super- 
fluid density. Here we consider the effective potential 
Uk as a function of the two independent variables n and 
/i. The flrst of equations P5)) states that in a Galilean 
invariant superfluid at zero temperature, the superfluid 
density is given by the full density of the fluid . Equa- 
tions (|43|) also imply that the Goldstone mode velocity Ck 
coincides with the macroscopic sound velocity [7,, J7,, 2^ > 
i.e. 



duk 
djjL 



nk 



(44) 



Since thermodynamic quantities, including the conden- 
sate "compressibility" dno.k/d^ should remain flnite in 
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Hpclk) 

FIG. 2: (Color online) A^, Zc,k and Va.u vs \n{j>G/k) where 
PG = a/ {gm)^n/4:TT for n = 0.01, 2mg — 0.1 and d — 2 
[ln(pG/pc) — —5.87]. The inset shows pa vs 2mg obtained 
from the criterion Va,pc = VA,k=o/2 [the Green solid line is 
a fit to PG ~ (2mg)^/^l 




lii{pG/k) 



not expected to flow in the infrared limit, they can be 
approximated by their fc = values. As for the longitu- 
dinal correlation function, its value is obtained from the 
replacement — ^ C(aj^ + c^p2)(3-'^)/2 (with C a con- 
stant). From (|¥T|l and we then deduce the exact 
infrared behavior of the normal and anomalous propaga- 
tors (at k = 0), 



Gnip) = - 



nomc 



1 



n 



2t-.2 



mc^ dno iu! 



(46) 



Ga„(p) 



Home 



n 



c2p2 



2Gii(p), 



where 



Gii(p) 



2noC(a;2 + c2p2)(3-d)/2- 



(47) 



The hydrodynamic approach of Sec. Illll correctlv predicts 
the leading terms of (|46p but approximates dno/dfi by 
n/mc^. On the other hand, it gives an explicit expres- 
sion of the coefficient C in the longitudinal correlation 
function (|T7)) . 



FIG. 3: (Color online) Condensate density no,fe, superfluid 
density ns,k and Goldstone mode velocity Ck vs \n{pG/k). The 
parameters are the same as in Fig. (2] 



B. RG flows 



the limit fe — )■ 0, we deduce from (|43p that Zc^fc ^ Afc '-^ 
vanishes in the infrared limit, and 

lim Cfc lim f ^^-^ ^ ^ . (45) 
fc^o \2mVA,k ) 

Both Zyi^fc = nt^/nQ_k and the macroscopic sound velocity 
Cfe being finite at fc = 0, VA,k (which vanishes in the 
Bogoliubov approximation) takes a non-zero value when 
fc -> 0. 

The suppression of Zc,k, together with a finite value 
of VA,k=o shows that the average effective action ([55)) 
exhibits a "relativistic" invariance in the infrared limit 
and therefore becomes equivalent to that of the classi- 
cal 0(2) model in dimensions d + 1. In the ordered 
phase, the coup ling constant of this model vanishes as 
Afe ^ ;j4-(d+i) jj^^ which is nothing but our starting 
assumption (HUl) . For fc — )■ 0, the existence of a linear 
spectrum is due to the relativistic form of the average 
effective action (rather than a non-zero value of Xkno^k 
as in the Bogoliubov approximation) [36j . 

To obtain the fc = limit of the propagators (at 
fixed p), one should in principle stop the flow when 
fc ~ a/p^ + up' I (? . Since thermodynamic quantities are 



The conclusions of the preceding section can be ob- 
tained more rigorously from the RG equation (|36p sat- 
isfied by the average effective action. The flow of Afc, 
Zc,k and VA,k is shown in Fig. [2] for a two-dimensional 
system in the weak-coupling limit. We clearly see that 
the Bogoliubov approximation breaks down at a charac- 
teristic momentum scale pQ ^ ^ {gm)-^n. In the Gold- 
stone regime k <^ pc, we find that both A^ and Zc,k 
vanish linearly with fc in agreement with the conclusions 
of Sec. IIV Al Furthermore, VA,k takes a finite value in 
the limit fc — > in agreement with the limiting value 
(|45l) of the Goldstone mode velocity. Figure [3] shows 
the behavior of the condensate density no,fe, the super- 
fiuid density Us^k — ■^A,fe?^o,fe and the velocity Cfe. Since 
ZA,k=o — 1.004, the mean boson density n^, = n^ fc 
is nearly equal to the condensate density rto,fc. Apart 
from a slight variation at the beginning of the flow, no.fc, 
ris.k — ZA.k^o^k and Ck do not change with fc. In par- 
ticular, they are not sensitive to the Ginzburg scale pc. 
This result is quite remarkable for the Goldstone mode 
velocity Ck, whose expression involves the param- 
eters Afc, Zc,k and VA,k, which all strongly vary when 
fc ~ PG- These flndings are a nice illustration of the fact 
that the divergence of the longitudinal susceptibility does 
not affect local gauge invariant quantities [l7l[20|. 
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0.01 0.02 0.03 0.04 0.05 0.06 0.07 
t/U 

FIG. 4: Phase diagram of the two-dimensional Bose-Hubbard 
model obtained from the NPRG [39] showing the first Mott 
lob corresponding to a density n — 1. 



V. THE SUPERFLUID-MOTT-INSULATOR 
TRANSITION 

The Bose-Hubbard model describes bosons on a d- 
dimensional lattice with Hamiltonian 13711 



H = - (V'i'0r' + h.c.) - ^ fir 

(r,r') 

— ^nr(nr - 1), 



where tj)'^'^ is a creation/annihilation operator defined at 
the lattice site r, (r, r') denotes nearest-neighbor sites, 
and rir = V'J^r- This model can be studied within the 
NPRG framework using a formulation adapted to lattice 
models (ssl . [39| . Here we focus on the critical behavior 
at the superfluid-Mott-insulator transition. In the low- 
energy limit, we can take the continuum limit and express 
the average effective action as 

TkW^]^ dr j d''r[<l)*{Zc,kdr-VA,kd'i 

- ZAMV'')4> + Vk{n)\, (49) 

where n(r,T) = |0(r, r)p and 0(r,T) = (■0(r, r)) denotes 
the superfluid order parameter. Near the transition, we 
can expand the effective potential Vk{n) about n — Q, 



(48) 



Vfe(n) 



14(0) + 4n + yn^. 



(50) 



The transition line in the {t/U, ^/U) plane is defined by 
5 = 6k=o = (Fig. 111). The condensate density no = 
no,fc=o is obtained from the minimum of V{n) = Vk=a{n). 
In the Mott phase (<5 > 0), tiq vanishes while no — —S/X 
in the superfluid phase ((5 < 0). The mean boson density 
(i.e. the mean number of bosons per site) is given by 



dV{0) 

dV{0) 
dfj. 



A dfi 



d\ 
2A^ dfj. 



(insulator), 
(superfluid) . 



(51) 





ZA.k 


VA.k 


Zc.k 




no.k 


superfluid 


z\ 


vx 


k 


k 


no 


critical 
behavior 


XY 






k 




k^+^ 


mean-field 


z\ 


vx 


Z*c 


|lnfc|~i 


P\\nk\-^ 


insulator 


z\ 


vx 


Z*c 


A* 






TABLE I: Critical behavior at the superfluid-Mott-insulator 
transition and infrared behavior in the superfluid and Mott- 
insulator phases (d = 2). The stared quantities indicate 
nonzero flxed-point values and rj denotes the anomalous di- 
mension at the three-dimensional XY critical point. 



In the Mott phase, the density —dV{0) /dfi is pinned to 
an integer value. 

The existence of two universality classes at the quan- 
tum phase transition between the superfluid and the 
Mott insulator follows from symmetry arguments (4dj . 
The last of the Ward identities (|3S|), which is associated 
to (local) gauge invariance [131 1 can be expressed as 



dS 
dfi 



S_dX 



(52) 



in the superfluid phase. At the tip of the Mott lob where 
both 6 and dS/d^ vanish, Zc vanishes. The effective 
action F^^o then exhibits a "relativistic" invariance and 
we expect the transition to be in the universality class 
of the {d + l)-dimensional XY model with a dynamical 
exponent z = 1. According to ([5T|l . the lob tip {S = 
dS/dfi = 0) corresponds to a transition taking place at 
constant density. Away from the tip, the transition is 
accompanied by a change in density. Zc is finite and the 
transition is mean-field like for d > 2 with z ~ 2 (with 
logarithmic corrections for d = 2) |40l] . 

Table U summarizes the results obtained for the two- 
dimensional Bose-Hubbard model [so']. The NPRG pro- 
vides a natural explanation for the critical behavior in the 
XY universality class which is observed when the tran- 
sition takes place at constant density (tip of the Mott 
lob). Indeed, since the infrared behavior in the super- 
fluid phase is characterized by a relativistic symmetry, 
it is not surprising that this symmetry remains at the 
transition to the Mott insulator. By contrast, the Bo- 
goliubov fixed point, which has a dynamical exponent 
z = 2, is clearly a poor starting point to understand the 
superfluid-Mott-insulator transition at the lob tip. 



VI. CONCLUSION 

Interacting bosons at zero temperature are character- 
ized by two momentum scales: the healing (or hydrody- 
namic) scale pc and the Ginzburg scale pc- For momenta 
IpI Pc, it is possible to use a hydrodynamic descrip- 
tion in terms of density and phase variables. This de- 
scription allows us to derive the correlation functions of 
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the order parameter field without encountering infrared 
divergences. In the Goldstone regime |p|, |w|/c pq, 
density fluctuations play no role any more and both the 
transverse and longitudinal correlation functions are fully 
determined by transverse (phase) fluctuations. In this 
momentum and frequency range, the coupling between 
transverse and longitudinal fluctuations leads to a diver- 
gence of the longitudinal susceptibility and singular self- 
energies. Nevertheless, in the weak-coupling limit, the 
Bogoliubov theory applies to a large part of the spec- 
trum where the dispersion is linear (|p| < p^) and breaks 
down only at very small momenta |p| pc Pc- More- 
over, thermodynamic quantities are not sensitive to the 
Ginzburg scale pc and can be deduced from the Bogoli- 
ubov approach. 



A direct computation of the order parameter correla- 
tion function (without relying on the hydrodynamic de- 
scription) is possible, but one then has to solve the prob- 
lem of infrared divergences which appear in perturbation 
theory when |p| < and signal the breakdown of the 
Bogoliubov approximation. The NPRG provides a nat- 
ural framework for such a calculation. It shows that in 
the Goldstone regime |p|, \uj\/c <^ pc, the system is de- 
scribed by an effective action with relativistic invariance 
similar to that of the {d + l)-dimensional classical 0(2) 
model. This similarity sheds light on the critical behavior 
of the superfluid-Mott-insulator transition in the Bose- 
Hubbard model which belongs to the XY universality 
class when the transition takes place at fixed density. 
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